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ABSTRACT 

The  objective  of  this  research  project  was  to  build  on  the 
earlier  work  of  BOURDEAU(1),  who  brought  together  the 

probabilistic  concept  of  stress  transmission  in  granular 

•  (2)  * 
bodies  (HARR  )  and  the  analogous  stochastically  based 

•  •  (3) 

deformation  mechanism  proposed  by  LITWINISZYN  ,  to 
investigate  the  response  of  particulate  material  to 
vertically  applied  slowly  moving  surface  loadings. 

A  thorough  review  of  existing  and  current  literature  on  the 
mechanics  of  particulate  media  yielded  the  discussions  on 
Porosity  (Chapter  Two) ,  Force  and  Stress  (Chapter  Three)  and 
Deformation  (Chapter  Four) .  As  HARR  and  LITWINISZYN 
utilized  an  explicitly  probabilistic  approach,  which  has 
been  adequately  confirmed  in  the  testing  and  coupled  model 
of  BOURDEAU,  it  was  decided  to  approach  the  current  project 
from  the  same  basis.  The  perceived  advantage  of  this 
methodology  is  that  the  response  of  granular  materials  may 
be  characterized  by  a  very  few  parameters,  whereas  the  older 
stress-dilatancy  and  yield  theories'4’ <5)  <6)  l7)  abound  in 
extra  parameters. 

The  testing  phase  was  undertaken  at  the  soil  mechanics 
laboratory  of  the  Swiss  Federal  Institute  of  Technology 
(EPFL)  at  Ecublens,  Switzerland.  A  fine  dry  silt  soil  was 
tested  under  both  repeated  static  loading  and  quasi-mobile 
loading.  An  array  of  lead  shot  was  placed  periodically 
throughout  the  soil  body,  and  records  were  made  of  the 
deformation  of  the  soil  under  the  surface  applied  loading 
through  exposure  of  the  soil  model  to  X-rays.  Load  was 
applied  to  the  surface  of  the  soil  through  a  rigid  rough 
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plate.  In  the  static  series  of  loads,  the  plate  was 
repeatedly  applied  to  the  same  point  on  the  soil  surface, 
whereas  in  the  moving  load  series,  the  load  plate  was 
relocated  contiguous  to  its  prior  position  at  each 
application.  In  this  way,  for  the  latter,  a  •’quasi-moving" 
load  mechanism  was  modeled.  Chapter  Five  details  the  test 
procedures  and  the  testing  sequence.  The  methodology  of 
data  reduction  is  also  outlined. 

Chapter  Six  provides  a  detailed  analysis  and  discussion  of 
the  results.  The  conclusion  is  drawn  that  the  response  of 
dry  granular  material  to  vertically  applied  moving  loads  is 
one-dimensional,  and  the  magnitude  of  deformation  may  be 
estimated  from  an  equivalent  static  loading.  There  is 
evidence  that  coefficient  of  lateral  stress  (v)  is  not 
constant  throughout  the  material  under  load  but  varies 
significantly  with  depth.  The  pattern  of  the  variation  of  v 
with  depth  under  the  load  is  not  inconsistent  with  HARR's 
modification  of  SCHMERTMANN ' s  method  of  settlement  analysis 
and  lends  credence  to  it. 

Of  all  the  data  analyzed,  it  is  evident  that  the  primary 
data,  that  is  the  horizontal  and  vertical  deformations,  are 
the  best  behaved.  The  magnitudes  of  the  observed  vertical 
deformations  are  generally  an  order  of  magnitude  greater 
than  the  observed  horizontal  deformations.  The  secondary, 
or  derived,  parameters  discusssed,  i.e.,  the  horizontal 
vertical  strains,  are  less  informative  and  somewhat  less 
well  behaved,  while  the  tertiary  parameters  (volumetric 
strain  and  porosity)  are  considered  only  qualitatively  since 
error  propagation  from  the  raw  deformation  data  may  be 
considerable. 
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Chapter  One 
INTRODUCTION 
Particulate  Hedja 

The  classic  continuum  of  mechanics  is  distinguished  by  the 
absence  of  voids  and  hence  complete  continuity  of  matter  is 
assumed  from  point  to  point  within  the  body  under 
consideration.  The  form,  or  shape  of  the  body  is  defined  by 
the  geometric  characteristics  of  its  boundary.  Inclusions 
or  voids  may  be  treated  explicitly  in  the  form  of  additional 
boundary  conditions,  with  the  necessary  a  priori  knowledge 
of  location,  shape,  and  size.  At  the  microscopic,  or 
molecular,  level,  continue  are  collections  of  discrete 
particles  (atoms  and  molecules) ,  wherein  the  inter-atomic 
and  inter-molecular  forces  are  sufficiently  strong  to  bind 
the  individual  elements  together  such  that  they  behave  as  a 
common  mass. 

Granular  or  particulate  materials,  on  the  other  hand,  are 
generally  characterized  by  being  collections  of  individual 
particles;  generally,  the  visible  scale.  Any 

inter-particle  attraction,  in  the  absence  of  moisture,  is 
sufficiently  weak  as  to  be  sensibly  neglected  in  engineering 
terms.  The  whole  mass  of  a  dry  .granular  material  is 
distinguished  from  a  continuum  by  being  "full  of  holes". 
The  individual  particles,  whether  spherical,  cubical  or 
irregular,  cannot  fit  together  perfectly,  but  rest  in 
point-to-point  contact  with  their  neighbors.  This  results  in 
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included  volumes  that  are  partially  empty,  or  void. 

Unlike  continua,  the  voids  in  granular  materials  cannot  be 
treated  explicitly:  they  are  so  random,  or  chaotic,  in 
location,  shape  and  size  as  to  defy  rational  description. 
Equally,  the  non-void,  or  solid  component  of  granular 
materials  is  not  spatially  continuous;  that  is,  in  spite  of 
the  existence  of  particle-to-particle  contact,  there  is  no 
continuity  in  all  directions  at  all  points. 

Force  and  Stress 

Classical  continuum  mechanics  uses  the  stress  at  a  point  as 
a  measure  of  the  transmission  of  boundary  energy  within  and 
through  the  continuum.  This  concept  arises  from  the 
infinitesimal  mathematics,  wherein  force  is  applied  over  an 
area,  which  becomes  stress  in  the  limit  as  the  area 
approaches  zero.  Consequently,  stress  must  be  defined  at  a 
point.  It  is  an  areal  parameter  taken  as  the  averaged  force 
in  the  vicinity  of  the  point.  In  continuum  mechanics,  this 
is  generally  brushed  aside  by  invoking  statistical 
macroscopic  equivalence;  however,  it  begins  to  take  on  very 
important  physical  meaning  when  applied  to  granular 
materials. 

BOUSSINESQ's  classical  solution  of  a  homogeneous, 
semi-infinite  elastic  continuum  has  often  been  applied  to 

.  ..  (Q) 

granular  materials.  FOPPL  and  FOPPL  ,  koegler  and 
SCHEIDIG(l0),  and  TURNBULL,  MAXWELL  and  AHLVIN(U)  among 
many  others  have  all  observed  and  commented  at  length  on  the 
inability  of  classical  mechanics  to  predict  adequately  the 
observed  response  of  granular  materials  to  loading. 

In  the  early  seventies,  also  concerned  by  the  shortcomings 
of  the  classical  solution,  the  principal  investigator 
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(HARR121)  developed  a  probabilistic  approach  based  upon  the 
stochastic  nature  of  the  propagation  of  force  into  a 
granular  material.  The  development  was  predicated  on  the 
application  of  the  methodology  of  a  random  walk. 


In  classical  continuum  mechanics,  the  deformation  of  a  body 
under  load  is  assumed  to  arise  from  the  induced  volume 
changes  and  distortion  of  the  solid  material :  there  are  no 
voids.  This  is  modeled  mathematically  by  ensuring  strain 
compatibility ,  or  continuity ,  and  by  invoking  a  constitutive 
law,  such  as  linear  elasticity  (Hooke's  Law) .  In 
particulate  media,  there  is  no  continuity  in  the  classical 
sense  as  the  particles  are  not  bound  to  each  other,  so  that 
one  particle  may  move  (translate,  rotate)  relative  to 
another  without  necessarily  affecting  its  neighbors. 
Further,  deformation  in  granular  materials  results,  not  from 
the  deformation  of  the  solid  material,  but  from  relative 
geometric  displacements  between  particles  ( BOURDEAU n 1 ) . 


’reject  Objective 


Building  upon  the  work  of  previous  studies,  the  objective  of 
chc  current  project  was  to  gain  a  better  understanding  of 
the  response  of  dry  cohesionless  granular  material  to  slowly 
moving  vertically  applied  surface  loadings. 


4 


Chapter  Two 

POROSITY 


The  classic  continuum  of  mechanics  is  distinguished  by  the 
absence  of  voids  and  hence  complete  continuity  of  matter  is 
assumed  from  point  to  point  within  the  body  under 
consideration.  The  form,  or  shape  of  the  body  is  defined  by 
the  geometric  characteristics  of  its  boundary.  Inclusions 
or  voids  may  be  treated  explicitly  in  the  form  of  additional 
boundary  conditions,  with  the  necessary  a  priori  knowledge 
of  location,  shape,  and  size.  At  the  microscopic,  or 
molecular,  level,  continua  are  collections  of  discrete 
particles  (atoms  and  molecules) ,  wherein  the  inter-atomic 
and  inter-molecular  forces  are  sufficiently  strong  to  bind 
the  individual  elements  together  such  that  they  behave  as  a 
common  mass. 

Granular  or  particulate  materials,  on  the  other  hand,  are 
generally  characterized  by  being  collections  of  individual 
particles;  generally,  at  the  visible  scale.  Any 
inter-particle  attraction,  in  the  absence  of  moisture,  is 
sufficiently  weak  as  to  be  sensibly  neglected  in  engineering 
terms.  The  whole  mass  of  a  dry  granular  material  is 
distinguished  from  a  continuum  by  being  ’’full  of  holes". 
The  individual  particles,  whether  spherical,  cubical  or 
irregular,  cannot  fit  together  perfectly,  but  rest  in 
point-to-point  contact  with  their  neighbobors.  This  results 
in  included  volumes  that  are  partially  empty,  or  void. 

Unlike  continua,  the  voids  in  granular  materials  cannot  be 


5 


treated  explicitly:  they  are  so  random,  or  chaotic,  in 
location,  shape  and  size  as  to  defy  rational  description. 
Ect  ally,  the  non-void,  or  solid  component  of  granular 
materials  is  not  spatially  continuous;  that  is,  in  spite  of 
the  existence  of  particle-to-particle  contact,  there  is  no 
continuity  in  all  directions  at  all  points. 


A  quantitative  measure  of  the  pore-structure  of  particulate 
materials  is  nrovided  by  their  "POROSITY",  which  is  defined 
classically  as  "  the  volume  of  voids  in  the  material,  as  a 
proportion  of  the  total  material  volume".  Thus: 


n 

V 


V  V 

V  V 


V  +  V  V 

V  B  T 


(1) 


where 


ny  is  the  porosity  (in  volume  terms) , 

V  is  the  volume  of  voids  in  the  material, 

V 

V  is  the  volume  of  solids  in  the  material,  and 

S 

V  is  the  total  volume  of  the  material  =  V  +  V 

T  vs 


Given  a  unit  cubic  volume  of  particulate  material,  each 
point  in  the  cube  may  be  assigned  coordinates,  x,  y,  z.  By 
sampling  in  the  volume  at  a  large  number  of  random  points 
IR(x, y,z),  and  noting  a  "hit"  each  time  that  point  is  in  a 
void  and  a  "miss"  in  the  case  of  a  solid,  it  will  be  noted 
that  the  ratio  -j— ■■  —  will  approach  a  constant. 
The  constant  value  about  which  the  sampling  ratio  fluctuates 
is  referred  to  as  the  "probability  of  a  hit"(42, .  This 
definition  of  porosity  may  be  interpreted  as;  porosity  is 
the  probability  of  finding  a  void  in  the  material,  or: 


n 

V 


#  hits 

#  hits  +  #  misses 


=  Probability 


“  - 
void 


(2) 


From  equation  (1),  it  will  be  seen  that,  mathematically, 

porosity  may  range  from  zero  (Vy  =  0  =  continuum)  to  unity 

fV  =  0  =  no  material!).  However,  normal  engineering  soils 
8 
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fall  into  a  narrower  range,  0.2  s  s  0.5,  depending  upon 
the  material  gradation  and  degree  of  compaction. 

(2)  . 

HARR  points  out  that  "  ideal"  arrangements  of  identical 
spheres  yield  the  following  theoretical  porosities: 


Packing  Arrangement  Theoretical  Porosity, 

Cubic  0.4764 

Orthorhombic  0.3954 


Tetragonal -spheroidal 
Rhombhedral 

TERZAGHI  and  PECK(36)provide 
porosities  for  real  soils: 

Soil  Description 
Loose  uniform  sand 
Dense  uniform  sand 
Loose  graded  sand 
Dense  graded  sand 
Glacial  Till 

SCHULTZEU0)  gives  a  similar 
including  their  standard  devial 
degree  of  compaction: 


0.3019 

0.2595 

the  following  typical 

Typical  Porosity, n^ 

0.46 

0.34 

0.40 

0.30 

0.20 

list  of  typical  values 
ons,  without  indication  of 


^oil  Description  Typical  Porosity ,ny  Standard  Deviation 


Gravelly  sand 

0.330 

0.062 

Coarse  sand 

0.399 

0.039 

"  dium  sand 

0.404 

0.041 

' xne  sand 

0.390 

0.032 

silt 

0.380 

0.051 

The  practical  upper  limit  for  porosity,  assuming  that  the 
particles  are  in  contact  with  at  least  two  others  (i.e.,  in 
their  least  stable  configuration)  is  about  0.60;  with 
arching  in  very  loose  materials,  this  might  reach  0.65.  The 
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lower  limit  of  porosity  is  less  well  defined,  and  depends 
largely  on  the  gradation  of  the  particles,  but  it  might  be 
taken  to  be  of  the  order  of  0.15. 

Porosity  is  not  a  material  parameter,  but  is  a  state 
parameter.  A  truck  load  of  granular  base  material  might 
have  a  porosity  of  0.40  when  freshly  tipped,  but  after 
compaction  it  may  well  have  a  porosity  of  0.30.  The 
material  has  not  changed  but  its  porosity  has! 

Given  that  a  granular  material  is  in  its  loosest  state,  it 
is  in  a  delicate  state  of  equilibrium.  The  only  way  that 
the  porosity  can  be  changed  is  by  applying  energy  to  the 
material  to  produce  a  re-arrangement  of  the  particles.  In 
general  terms,  a  granular  material  will  achieve  a  denser, 
more  stable,  configuration  (if  only  locally)  upon  the 
application  of  energy  to  it,  in  general  compliance  with  the 
principle  of  maximal  entropy  (2nd  law  of  thermodynamics) . 

(41 ' 

GRIVAS  and  HARR  point  out  that  whatever  the  packing 
arrangement  of  the  soil  particles,  the  general  measure  of 
particle  stability  is  given  by  the  coordination  number,  N, 
which  is  the  number  of  particle-to-particle  contacts  per 
particle.  The  greater  the  coordination  number,  the  more 
stable  the  packing,  in  its  state  of  local  equilibrium.  They 
point  out  that  porosity  and  coordination  number  appear  to  be 
closely  related  by  the  relationship  n/N.  »  n.  (The  value  n 
is  not  a  computed  result,  but  the  observed  value  was 
sufficiently  close  and  the  use  of  n  lends  a  certain  degree 
of  technical  elegance) . 

"Porosity  at  a  point"  has  no  meaning.  At  a  point  there  is 
either  a  void  ...or  a  solid.  Thus  porosity  is  a  statistical 
property  of  the  material,  and  if  it  is  used  in  reference  to 
a  point  location  within  the  material,  it  must  be  understood 
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that  it  refers  to  a  sufficient  volume  of  the  material 
surrounding  the  point  to  render  the  measure  s cable. 

Porosity  is  relatively  easy  to  measure  in  the  laboratory,  at 
least  in  a  bulk  sense,  once  it  is  assumed  to  be  constant 
throughout  the  sample.  In  the  field,  porosity  may  be 
estimated  from  "density  holes",  but  once  again,  it  must  be 
assumed  that  it  is  uniform  throughout  the  sample. 

In  the  course  of  the  present  research  project,  X-ray 
techniques  were  used  ( vide  Chapter  5)  which  permitted  an 
estimation  of  the  porosity  within  discrete  areas  of  the  test 
material.  As  will  be  shown  subsequently,  it  was  assumed 
that  the  initial  porosity  was  uniform  throughout  the 
material,  i.e.,  that  the  material  was  statistically 
homogeneous. 

Other  descriptors  exist  that  scale  the  state  of  voids  within 
granular  media.  The  void  ratio,  e,  which  is  the  ratio  of 
the  total  volume  of  voids  to  the  total  solid  volume  (a 
transformation  of  the  porosity).  The  relative  density,  Dr, 
purports  to  give  a  measure  of  the  state  of  the  'denseness' 
of  the  medium,  but  it  relies  on  attaining  measures  of  the 
loosest  and  densest  densities  achievable  in  empirical 
laboratory  tests. 

The  concept  of  porosity  is  central  to  any  competent 
description  of  granular  materials.  There  is  no  other 
parameter  which  so  well  describes  the  state  of  a  soil  and 
yet  is  fundamental.  The  understanding  of  porosity  and  of 
its  changes  is  the  sine  qua  non  of  the  mechanics  of 
particulate  materials. 
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Chapter  Three 
FORCE  and  STRESS 

Classical  continuum  mechanics  uses  the  stress  at  a  point  as 
a  measure  of  the  transmission  of  boundary  energy  within  and 
through  the  continuum.  This  concept  arises  from  the 
infinitesimal  mathematics,  wherein  force  is  applied  over  an 
area,  which  becomes  stress  in  the  limit  as  the  area 
approaches  zero,  consequently,  stress  is  defined  at  a  point, 
it  is  an  areal  parameter  taken  as  the  averaged  force  in  the 
vicinity  of  the  point.  In  continuum  mechanics,  this  is 
generally  brushed  aside  by  invoking  statistical  macroscopic 
equivalence;  however,  it  begins  to  take  on  very  important 
physical  meaning  when  applied  to  granular  materials. 

BOUSSINESQ ' s  classical  solution  of  a  homogeneous, 

semi-infinite  elastic  continuum  has  often  been  applied  to 

(9) 

granular  materials.  FOPPL  and  FOPPL  ,  KOEGLER  and 
SCHEIDIG(11),  and  TURNBULL,  MAXWELL  and  AHLVIN(12)  among 
many  others  have  all  observed  and  commented  at  length  on  the 
inability  of  classical  mechanics  to  predict  adequately  the 
observed  response  of  granular  materials  to  loading. 
Marginal  improvements  have  been  gained  by  assuming 

/  Ip) 

anisotropic  properties  (GARRARD  and  MULHOLLAND  ,  VAN 
CAUWELAERT(13) )  ,  and  by  providing  empirically  derived 
stress-dependent  moduli  (DUNCAN  and  CHANG<14),  JANBU(15), 
and  LADE065). 

In  1937,  POKROVSKI <17)  obtained  a  gaussian  type  of 
distribution  for  stresses  in  a  granular  soil  from  a  purely 
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deterministic  approach.  Following  this  work,  and  perhaps 
prompted  by  it,  KOGLER  and  SCHEIDIG<10)  suggested  that  a 

probabilistic  approach  might  be  justified.  The 

(18) 

probabilistic  approach  was  first  taken  up  by  KANDAUROV 
and  MULLER  ,  although  they  treated  the  results 
deterministically;  not  admitting,  or  recognizing  the 
inherent  stochastics  of  the  system.  SERGEEV(20) ,  in  the 

late  60/s  was  the  first  to  apply  direct  probabilistic 
reasoning  based  upon  Markovian  processes.  Independently, 
and  at  about  the  same  time,  SMOLTCZYK  was  achieving  the 
same  results,  but  with  a  much  less  rigorous  approach. 

In  the  early  seventies,  HARR  developed  a  probabilistic 
approach  based  upon  the  stochastic  nature  of  the  propagation 
of  force  into  a  granular  material,  which  also  introduced  the 
concept  of  the  distribution  of  stress  at  a  point.  His 

development  was  predicated  on  the  application  of  the 

(22) 

methodology  of  a  random  walk.  ENDLEY  and  PEYROT  ,  about 
the  same  time,  presented  the  concept  based  on  the  stiffness 
of  stress  paths.  In  1984,  GOLDEN  published  a  thorough 
examination  of  the  stochastic  theory  of  stress  propagation. 
While  he  did  not  add  materially  to  the  work  of  HARR,  his 
more  rigorous  development  and  critical  examination  of  the 
basic  assumptions  provided  a  measure  of  conf irmation  and 
support  to  the  general  probabilistic  approach. 

It  should  be  noted  that  the  brief  historical  outline 
presented  above  deals  only  with  the  concept  of  stress.  No 
implication  or  consideration  of  deformation  from  a  purely 
probabilistic  basis  was  discussed,  or  detailed,  by  the  noted 
authors.  When  deformation  was  considered,  it  followed  from 
the  introduction  of  a  linearization  of  Hooke's  law.  Thus 
the  probabilistic  theory  of  stress  propagation  in  granular 
materials  has  been  developed  in  a  vacuum  (dare  one  say 
it . . . in  a  void) . 
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Contact  Forces  in  Granular  Soils 

As  was  pointed  out  above,  the  deterministic  concept  of 
stress  is  a  convenient  artifice.  The  probabilistic  concept  will  be 
developed  below  from  the  starting  point  of  force,  which  is  a 
fundamental  physical  reality  and  is  readily  measured. 

The  subsequent  development  borrows  essentially  from  the 
earlier  work  of  HARR  .  As  with  the  works  of 

KANDAUROv‘18>  ,  HULLER<19>  and  GOLDEN<23),  it  assumes  that 
the  granular  medium  is  statistically  homogeneous  and 
isotropic  and  remains  so  throughout  the  loading  process. 

If  a  unit  point  load  is  applied  vertically  to  the  horizontal 
surface  of  a  granular  medium,  then  purely  from  consideration 
of  static  equilibrium,  the  applied  load  (force)  is  divided 
among  the  particles  upon  which  the  loaded  particle  rests. 

This  process  is  repeated  by  those  supporting  particles  to 
the  next  "layer”  of  particles.  The  precise  determination  of 
these  reactions,  their  magnitudes  and  directions,  is 
impossible.  However,  given  the  homogeneity  of  the  material, 
and  invoking  the  principle  of  maximum  entropy,  the  least 
biased  division  is  50:50.  This  may  be  represented 
schematically  in  two  dimensions  as  shown  in  Figures  1  and  2. 

Mathematically,  these  (Figures  1  and  2)  may  be  represented 
by  the  recursion  equation: 

Fa(x,2+Az)  =  p*Fz(x-Ax,z)  +  q*F  z(x+Ax,z)  (3) 

where  p  +  q  =  1,  and  p,  q  represent  the  probabilities  of  the 
force  dividing  to  the  "left"  and  to  the  "right" 
respectively,  and  Fz(*)  denotes  the  contact  force  in  the 
z-direction,  acting  on  the  particle  at  the  noted  coordinate 
point  ( • ) 
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Figure 


F  (x-Ax,z) 
2 


F  (x+Ax,z) 


Figure  1  Resolution  of  particle  forces  in 
a  granular  medium 


•  1  ■  — * —  — 
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2  *^°?«bilis^iC  Pr°P?9ation  of  particle  force  in 
an  isotropic  granular  medium 
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When  p  =  q  =  1/2,  then  equation  (3)  reduces  to: 

F  (x,z+Az)  =  i*F  (x-Ax, z)  +  i*F  (x+Ax,z)  (4) 

z  2  z  2  z 


As  noted  equations  (3)  and  (4)  are  recurrence  equations,  the 
second  of  which  may  be  transformed  into  the  partial 
differential  equation: 


8Fz(x,z) 

dz 


1  8m 
2 ‘ax2 


D  *F  (X,z) 


(5) 


where : 


ID 


=  Lim 
Ax->o 
Az-*o 


This  is  the  classical  Fokker -Planck  equation. 


Further,  it  may  be  stated  that,  if  p  »  q,  then  the  material 
is  isotropic,  but  if  p  *  q,  the  material  exhibits  a  degree 
of  anisotropy.  GOLDEN  showed  that  the  case  where  p  *  q  may 
also  be  evidenced  in  an  isotropic  material  when  the  imposed 
force  has  a  horizontal  component. 


HARR*2’,  GOLDEN,23)  and  BOURDEAU*1’  further  assumed  that  the 
diffusion  parameter,  ID  ,  is  a  function  of  z  only,  which  then 
removes  it  from  the  brackets  under  the  differential  with 
respect  to  x.  In  this  case: 


SFz(x,z) 

5Z 


where  D*  is  a  function  of  z  only, 


or  D *(z) . 


(6) 


Equation  (6)  forms  the  basis  of  the  works  of  KANDAUROV(18) , 
MULLER(19>,  HARR<2),  BOURDEAU'11  and  GOLDEN* 23 ’ ;  and  has  the 
solution: 
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FJx.z)  =  —  -L-—  •  e  2D  z  (7) 

i/2TTID*2 

The  notation  F *(x,z)  denotes  the  expectation  of  the  force  at 
that  point. 

The  physical  meaning  of  ID*  may  be  understood  from  the 
observation  that  D  z  is  the  equivalent  of  the  variance  of 
the  gaussian  distribution  represented  by  Eq.  7,  and  as  such 
is  a  measure  of  the  spread  of  the  distribution  of  force  on  a 
horizontal  plane.  It  is  a  property  of  the  material  and  must 
necessarily  be  positive. 

In  order  that  tensile  forces  do  not  occur  (a  basic 
assumption  in  dry  particulate  materials)  and  also  that  the 
force  component  F^(x,z)  tends  to  zero  at  large  z,  it  is 
necessary  that  the  product  ID*z  approaches  zero  as  z 
increases.  GOLDEN<23>  shows  that  ID*  is  a  power  function  of 
z,  and  stated  that  without  further  knowledge,  the  linear 

function  (i.e.,  Harr's)  would  lead  to  a  dimensionally 

.  ,  .  (2) 

balanced  exponent  in  equation  7.  As  pointed  out  by  HARR 

(23)  ,  •  . 

and  GOLDEN  ,  the  simplest  choice  consistent  with  these 
requirements  is  to  make: 


D  =  VZ  (7) 


for  all  z,  where  v  is  a  constant,  which  Harr  called  the 
11  coefficient  of  lateral  stress" .  Note:  while  Harr 
originally  developed  this  concept  from  a  consideration  of 
stress,  the  development  here  is  slightly  different,  based  as 
it  is  on  force  rather  than  stress. 


HARR(2)  showed  that  the  expected  horizontal  component  of  the 
force  is: 
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a2?  ( x,z ) 

F  <x,z)  =  VF  (x,z)  +  v2z2 - i— -  (9) 

z  dx2 

Assuming  that  the  second  term  on  the  right-hand  side  is 
small  in  comparison  with  the  first  term;  say,  for  example, 
if  Fz(x,z)  is  a  linear  function  of  x,  then: 


F  (*,z) 

v  =  -2 - 

F2(x,z) 


(10) 


He  also  pointed  out  the  analogy  between  this  definition  of 
v  and  that  of  the  coefficient  of  lateral  earth  pressure  of 
classical  soil  mechanics,  K: 


K  = 


(11) 


For  natural  soils  under  homogeneous  overburden  conditions, 
the  coefficient  of  lateral  stress,  v,  may  then  be 
approximated  by  Kq,  the  coefficient  of  earth  pressure  at 
rest,  such  that: 


D*  =  K0-Z  (12) 

Using  the  same  analogy,  HARR  proposed  that  u  be  bounded 
by  active  and  passive  analogues.  He  proposed  using  the 
REIMBERT  and  REIMBERT(24)  inequality: 


45  -  (3/2 

2 

£  £ 

f  ' 

45  +  E/2 

45  +  (3/2 

L  - 

45  -  E/2 

> 

where  )3  is  the  angle  of  repose  of  the  material  in  their 
formulation. 
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/pi 

HARR  and  GOLDEN  ,  when  discussing  v,  both  state  "..in 
a  homogeneous  medium..'',  while  BOURDEAUm  states  This 

important  simplification  is  the  same  as  neglecting  the 
effects  of  structural  modifications  induced  by  loading  on 
the  way  that  stresses  are  distributed  through  the  body  at 
equilibrium"  ( Author's  translation) .  While  this  subject  is 
the  topic  of  a  later  chapter,  the  basis  of  this  remark  stems 
from  the  belief  that  local  volumetric  deformations  are 
directly  related  to  changes  in  the  local  stress  field.  As 
the  force  field  represented  by  Eq.  (7)  is  obviously 
non-uniform  in  its  domain,  so  also,  the  porosity  of  the  body 
after  loading  would  no  longer  be  expected  to  be 
statistically  homogeneous.  It  is  advisable  to  state  its 
implications  at  this  point.  Making  the  assumption  that 
local  volume  changes  within  the  body  occur  as  the  result  of 
changes  in  the  local  forces,  and  that  deformation  of  the 
individual  particles  is  very  small,  if  not  negligible,  and 
therefore  that  deformation  proceeds  as  a  result  of  changes 
in  porosity,  (i.e.,  local  changes  in  the  geometric 
arrangements  of  the  particles)  ;  the  state  of  stress  at  a 
point  in  the  body  should  depends  upon  the  state  of  the 
material,  that  is,  the  local  stress  in  a  dense  material  (low 
porosity)  will  be  different  from  that  in  a  loose  material 
(high  porosity),  all  other  things  being  equal. 

In  the  previous  paragraphs  it  was  stated  that  deformation 

proceeds  from  changes  in  stress,  and  yet,  stress  has  not  yet 

(2) 

been  defined  within  the  present  context.  HARR  and 

(23)  •  •  • 

GOLDEN  both  treated  Eq.  7  as  a  stress  relationship;  in 
fact  HARR  used  the  notation,  s^,  for  the  expected  value 
of  stress,  so  that: 


sx(x,  z) 


v2nD*i 


2  D*z 


(14) 
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Another  derivation  of  stress  distribution  will  be  presented 
that  takes  the  material  state  into  account  more  explicitly. 

Stress 

Consider  a  small  volume  of  a  granular  medium  surrounding  the 
point  (x,z) ,  large  enough  to  define  a  competent  statistical 
sample  of  particles.  Within  that  sample  there  are  voids, 
which  cannot  sustain  force,  and  solids  which  accomodate  to 
all  induced  forces.  The  forces  acting  on  each  particle  are 
not  equal,  but  are  distributed  about  a  vertical  component  of 
a  mean  value,  Y^(x,z).  Thus,  we  may  write: 

Pro,  [void]  .  [f 0Vr;f  irioid]  +  **>  H  '  [forcf  Article] 

[vertical  1 

forces  averaged,  throughout  the  volume] 

Similar  relationships  may  be  written  for  forces  in  any 
geometric  orientation. 

As  was  shown  in  Chapter  Two,  the  probability  of  finding  a 
void  in  a  particulate  material  is  the  porosity,  nv.  Thus, 
as  there  are  no  forces  in  the  voids  of  a  dry  material: 

nv-0  +  (1  -  nv)-F z(x,z)  =  Sz(x,z)  (15) 

where,  S|x,z)  is  the  expected  value  of  the  expectation  of 
the  vertical  forces  in  the  volume  surrounding  the  point 

(x,z) • 

Now  from  equation  15,  it  may  be  seen  that: 

Sz(x,z)  =  (1  -  ny) -F ;z(x,z) 


(16) 
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which  shows  that  S^(x,z)  is  a  unique  value  for  each 
distribution  of  force  and  the  inherent  porosity  of  the 
material . 

It  is  proposed  that  the  variable  S z(x,z)  in  equation  16  be 
called  the  expected  value  of  local  stress  at  the  point 
( x,z )  .  It  is  seen  then  that  it  is  a  function  of  the  local 
material  state  (porosity),  whereas,  local  forces  are  not.  A 
material  body  under  loading  becomes  statistically 
nonhomogeneous ;  but  while  local  contact  forces  may  not 
change,  the  state  of  stress  will  reflect  changes  in 
porosity.  Further,  during  settlement  (or  consolidation) 
under  a  constant  applied  force,  the  local  stresses  and 
porosities  must  exhibit  the  relationship: 

S  (x,z) 

?  h~  )  ~  ConstaTVt  (17) 

'  V  ' 

So,  for  example,  if  a  material  exhibited  a  local  stress  of 
S  (x,z)  at  a  porosity  of  0.4,  and  after  some  time  the 
porosity  had  changed  to  0.35  under  the  same  loading 
conditions  (force  field),  the  local  stress  would  change  to 
1 . 083’S^,  an  increase  of  8%;  again  without  any  increase  in 
local  forces.  This  duality  of  local  porosity  and  local 
stress  is  one  of  the  central  ideas  of  the  present 
study . 
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Chapter  Four 
DEFORMATION 

In  classical  continuum  mechanics,  the  deformation  of  a  body 
under  load  is  assumed  to  arise  from  the  induced  volume 
changes  and  distortion  of  the  solid  material :  there  are  no 
voids .  This  is  modeled  mathematically  by  ensuring  strain 
compatibility,  or  continuity ,  and  by  invoking  a  constitutive 
law,  such  as  linear  elasticity  (Hooke's  Law) .  In 
particulate  media,  there  is  no  continuity  in  the  classical 
sense  as  the  particles  are  not  bound  to  each  other,  so  that 
one  particle  may  move  (translate,  rotate)  relative  to 
another  without  necessarily  affecting  its  neighbors. 
Further,  deformation  in  granular  materials  results,  not  from 
the  deformation  of  the  solid  material,  but  from  relative 
geometric  displacements  between  particles  ( BOURDEAU ( 1 1 , 
MARSAL125’)  . 

The  behavioral  differences  between  particulate  and  continua 
media  is  fundamental.  In  the  latter  case,  the  deformation 
of  an  elementary  volume  (3D)  or  area  (2D)  must  be  compatible 
with  those  elements  immediately  surrounding  it,  regardless 
of  the  adopted  constitutive  model.  This  is  a  result  of  the 
requirements  of  continuity,  or  compatibility.  In  elastic 
models  (favored  for  superposition  effects) ,  deformation 
follows  immediately  upon  application  of  load,  and  is 
immediately  reversible  upon  the  removal  of  that  load.  In 
other  models  (plastic,  visco-elastic  etc.),  this  is  not  a 
requirement,  but  then  the  benefits  of  superposition  are 
lost,  and  in  general,  there  is  an  increase  in  the  number  of 
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parameters  required  to  characterize  the  material  behavior. 

The  stochastically-based  models  of  deformation  in 
particulate  media  proposed  by  MARSALl2S> ,  LITWINISZYN<31  and 
BOURDEAU ' 1 }  do  not  require  'continuity ' ;  they  intrinsically 
incorporate  transience,  or  'time  effects'  and  are 
parsimonious  in  terms  of  parameters.  It  may  be  that  these 
models  are  overly  simplistic,  but  their  proven  advantage 
lies  in  the  observed  immediate  improvement  in  predictions  of 
steady-state  deformations. 

LITWINISZYN  (the  doyen  of  the  Polish  school)  developed  a 
comprehensive  stochastic  model  for  the  deformation  of 
particulate  media  resulting  from  volume  changes  within  the 
media  (e.g.,  the  deformation  resulting  from  the  collapse  of 
mining  galleries) :  upon  the  collapse  of  a  gallery,  the 
surrounding  rock  "falls"  into  the  void,  thereby  exchanging 
position  with  an  equivalent  volume  of  "space",  or  "void". 
It  is  this  progressive  mechanism  which  he  hypothesized  that 
causes  the  ultimate  surface  deformation  and  accounts  for  the 
transient  effect.  BOGDANOFF  and  SWEET(26)  extended  this 
work.  MARSAL<25>  made  the  connection  with  stressed  media 
(rather  than  'collapsed'  material)  ,  where  the  initiating 
internal  volume  change  resulted  from  local  changes  in  the 
stress  field.  B0URDEAU(1>  showed  that  while  MARSAL's<25) 
stochastic  model  was  organically  correct,  the  basic 
assumption  of  particle  'migration'  could  not  be  supported: 
he  proposed  the  coupling  of  HARR's  stress  propagation  model 
with  LITWINISZYN 'st3>  void  'migration'  deformation  model. 
BOURDEAU ( 1  *  treated  the  case  of  static  surface  loads  applied 
to  a  portion  of  the  surface  of  a  soil  body:  transience  was 
considered  only  as  component  of  the  steady-state 
deformation. 
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Deformation  of  Granular  Media 

There  are  three  primary  mechanisms  of  deformation  in  granular 
media <25> : 

i)  Elastic  deformation  of  the  solid  particles, 

ii)  Particle  crushing,  and 

iii)  Particle  realignment  and  relocation. 

Elastic  Deformation 

The  elastic  modulus  (E)  and  Poisson's  ratio  (u)  of  a  silica 
sand  sphere  (grain)  are  typically  6,000,000  lb/in2  and  0.15 
respectively  .  Uniaxial  deformation  under  a  unit 

diametrically  applied  load  will  be  of  the  order  of  0.0001 
inches  in  a  particle  0.04  inches  in  diameter,  or  an  overall 
strain  of  0.35%  .  Observed  deformations  of  assemblages 

of  granular  materials  generally  exceed  these  values  by 
orders  of  magnitude.  Undoubtedly,  in  a  loaded  granular 
medium,  there  will  be  a  component  of  elastic  particle 
deformation,  which  will  be  recoverable  upon  removal  of  the 
load.  The  magnitude  of  this  component  has  been  found  to  be 
negligible. 


Particle  Crushing 

The  expected  values  of  the  crushing  strengths  of  typical 
granular  materials  (in  civil  engineering  applications)  are 
of  the  order  of  30,000  lb/in2  (2160  ton/ft2) <27> .  Once 
again,  this  is  far  in  excess  of  the  expected  value  of 
stresses  induced  within  a  loaded  granular  mass.  While 
expected  values  of  stresses  in  granular  bodies  may  be 
relatively  low,  point  contacts  between  particles  can 
engender  extremely  high  local  stresses  sufficient  to  cause 
crushing.  This  mechanism  is  progressive  until  the  contacts 
have  degraded  to  such  an  extent  that  the  area  of  contact  has 
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increased  sufficiently  to  reduce  the  local  stresses  below 
critical  levels.  In  some  rock  and  granular  fills,  where 
the  quality  of  the  stone  is  either  variable  or  poorly 
controlled,  there  may  well  be  some  individual  particles 
which  undergo  crushing,  however,  in  general  practice,  this 
should  be  relatively  rare  .  While  there  will  be  some 
deformation  resulting  from  particle  crushing,  the  incidence 
of  significant  crushing  is  negligible. 

Particle  realignment  and  relocation. 

The  deformation  under  load  of  particulate  media  is  based 
upon  the  premise  that  in  a  state  of  equilibrium,  or 
steady-state,  the  forces  acting  on  each  particle  are  in 
balance  .  With  increases  in  the  imposed  loads,  local 

changes  occur  in  both  the  magnitude  and  direction  of  the 
forces  acting  on  the  particles.  The  departure  from  static 
equilibrium,  may  render  particles  unstable,  in  which  case 
they  may  move  to  a  more  stable  configurations.  If  the 
particles  in  a  local  volume  remain  in  undisturbed 
equilibrium,  there  is  no  change  in  the  volumetric 
distribution  of  solids  and  voids  in  that  region.  However, 
if,  in  order  to  acquire  a  state  of  enhanced  equilibrium,  the 
particle  slips  or  rotates,  then  there  may  result  a  local 
re-arrangement  of  the  particles.  As  BOURDEAUn>  and  others 
have  shown,  this  leads  to  a  redistribution  of  the  local 
solid  density.  If  the  motion  leads  to  a  consequent  local 
densif ication,  the  "excess"  voids  successively  exchange 
positions  with  those  of  particles  located  above.  That  is  to 
say,  the  "voids"  travel  upward  to  the  surface  in  much  the 
same  way  as  gas  bubbles  in  a  liquid. 

BOURDEAU(1)  pointed  out  that  if  a  single  particle  is 
removed  from  a  stable  particulate  arrangement  the  newly 
created  void  will  (can)  be  replaced  by  a  solid  particle 
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(thus  increasing  the  local  density) ;  but,  thereafter  on  its 
random  flight  through  the  material,  the  exchanged  "bubble" 
of  air  has  no  further  densifying  effect,  since  the  'before' 
and  'after'  situations  are  essentially  undisturbed. 
Figure  3  schematically  demonstrates  this  point.  Due  to  a 
change  in  local  stresses,  the  void  "particle"  a,  is  replaced 
by  solid  particle  b,  thus  increasing  the  local  density  and 
stability.  The  new  void  b*,  is  in  turn  displaced  by  solid 
particle  c,  and  so  on.  After  all  movement  has  ceased,  the 
local  density  in  frame  A  has  increased  due  to  the  capture  of 
particle  b  to  a,  while  the  local  density  in  frame  B  remains 
unchanged,  since  its  mass  balance  remains  constant,  i.e.,  it 
has  both  gained  a  particle  (e)  and  lost  a  particle  (c) . 


The  following  is  based  on  the  work  of  LITWINISZYN(3> ,  and 
provides  a  steady-state  model  of  the  effect  of  void 
migration.  Denoting  the  volume  of  a  void  displaced  at  point 
(x,z)  within  the  body  as  wq,  the  recurrence  relationship  for 
the  passage  of  the  void  to  the  surface  may  be  given  as: 


wq(x,z+Az)  =  §•  [w0(x-Ax,z)+wo(x+Ax,z)  j 
wo(x,z+Az)  -  wq(x,z)  =  \'  [wo(x-Ax,z)-2wo(x,z)+wo(x+Ax)  J 


(18) 

(19) 


awo(x,z)  x  9Wo(X,Z) 
—  ...  -  — — — —  =  -»Q»„  ■  —  - 

3z  2  „  _  _  2 


(20) 


where 


0  =  Lim 


3x 


t  2^ 
Ax 


Az 

V  / 


Again,  as  in  equation  (14)  the  solution  may  be  written  as: 


,  %  wo  "  2Qz 

wQ(x,z)  - - - 


v4tiQ  z 


(21) 
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This  relationship  provides  a  measure  of  the  location  of  the 
void  in  the  body  but  not  its  size.  The  initiating  void 
volume,  wq,  is  a  function  of  particle  sizes,  and  the 
skeletal  arrangement  of  the  grains  at  the  point  of 
origination  of  the  void.  To  make  the  orientation  of  the  "z" 
variable  in  the  stress  and  void  propagation  relationships 
compatible,  the  above  equation  may  be  rewritten: 


wQ(x,z) 


w. 


/27in  (h-z) 


X 

2fJ(h-z) 


(22) 


where  h  is  the  depth  below  the  surface  of  the  originating 
solid-void  exchange. 


BOURDEAU ( 1 1  defines  the  change  in  local  volume  from  a  local 
area  within  the  body  as: 


AV  -  ~*V 
P 


(23) 


where  ha 

AV 

V 

/3 


is  the  local  change  in  the  first  stress 
invariant, 

is  the  local  volume  change, 
is  the  initial  volume, 
is  a  "compression  modulus” 


In  this  case,  it  will  be  seen  that  the  local  volume  change 
AV,  is  the  sum  of  the  individual  voids  lost  to  the  volume, 
or: 

AV  =  £w‘(x,z) 

The  local  volumetric  strain  may  be  written  as: 


AV/V  *  ci  +  e2  +  c3 

where  the  ej  are  mutually  orthogonal  normal  strains. 


(24) 
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Pleading  linear  elasticity  with  respect  to  equation  (24),  it 
can  be  demonstrated  that  the  parameter  /3  in  equation  (23) 
may  be  expressed  in  terms  of  linear  elastic  constants,  E  and 
*i : 


For  2-D  (plane  strain) : 


c  +  c 

X  z 


+  a 

i 


(25a) 


or 


c 

X 


e 

z 


(1-2*0*  (ffx 


+  ffz> 


+  tr  ) 


thus 

B  _  E 

P  (1+U)  (l-2*i) 


(25b) 


(25c) 


For  3-D 


c  +  c  +  c  —  =•  (cr  +  0-  +  <r 

x  y  z  £  x  y  z 


2*1  (cr  +  cr  +  <j  )  ) 

x  y  z  ' 


(26a) 


or 


c  +  c  +  c 

x  y  z 


1-2*1 

E 


+  <T 

y 


CTz> 


+  a 

n 


+  cr 


iii 


) 

(26b) 


thus 


P  - 


E 

(1-2*1) 


(26c) 


It  is  arguable  that  these  relationships  hold  true  for  a 
granular  material,  but  in  any  case,  they  do  strike  an 
interesting  analogy. 


B0URDEAUa) ,  in  using  his  relationship,  implicitly  accepted 
deformation  as  a  continuous  reversible  and  linear  function 
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of  stress.  However,  upon  the  removal  of  load  from  such  a 
body,  deformation  is  not  fully  recovered,  (consider,  for 
example,  the  permanency  of  footsteps  in  sand) .  Upon  the 
replacement  of  the  load,  deform ition  is  essentially  linear 
up  to  the  previous  maximum  magnitude  of  load;  thereafter, 
deformation  is  non-linear  and  transient.  (Analogous  to 
precompression  and  virgin  compression  in  the  consolidation 
of  clays) .  Although  the  deformation  is  linear  in  the 
"elastic”  or  pre-loaded  range  when  the  surface  is  re-loaded, 
this  holds  true  only  when  the  load  is  re-applied  at  its 
original  location  on  the  surface.  When  the  load  is 
re-applied  at  a  different  surface  location,  the  induced 
stress  state  at  a  point  in  the  body  may  indeed  not  exceed 
the  previous  magnitude  of  the  trace  (o^  +  <r  +  <r  )  ;  but 
the  principal  directions  may  well  have  changed  orientation. 
This  possibility  is  not  accounted  for  in  equation  (23) . 

The  exhibited  "elastic"  portion  of  the  load  -  reload 
response  may  be  explained  as  follows.  Consider  the  state  of 
a  statistically  valid  volume  around  some  point  (x,z)  in  a 
granular  body  under  a  given  steady-state  stress  regime 
wherein  the  particles  are  all  in  static  equilibrium.  This 
equilibrium  is  not  necessarily  critical.  If  the  stress 
state  changes,  the  equilibrium  between  the  particles  may  be 
disturbed  and  a  permanent,  or  irreversible  deformation  will 
occur.  However,  it  is  wholly  possible  that  the  equilibrium 
state  prior  to  loading  is  sufficiently  robust  under  the 
increment  of  load  so  that  the  particles  retain  their 
relative  positions,  in  which  case  the  applied  energy  is 
absorbed  by  elastic  deformation  of  the  particles. 

Unlike  the  concept  of  stress,  deformations  are  directly 
measurable.  The  experimental  phase  of  the  present  project 
wilt  make  specific  use  of  this  advantage  to  monitor  the 
internal  deformations  of  granular  bodies  during  loading. 


28 


(2) 

In  summary,  HARR  developed  the  concept  of  stress  in 
granular  materials  without  consideration  of  deformation  and 

(3) 

LITWINISZYN  theorized  the  mechanism  of  deformation 
without  consideration  of  stress.  BOURDEAU(1>  was  the  first 
to  bring  together  the  two  disparate  hypotheses  within  a 
probabilistic  (rather  than  from  a  deterministic)  framework. 
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Chapter  Five 
EXPERIMENT 

The  experimental  phase  of  this  research  project  fell  into 
two  clearly  definable  parts:  the  response  to  repeated  static 
loadings  and  to  incremental  moving  loadings.  Consequently, 
the  initial  experiments  were  set  up  to  identify  and  quantify 
material  parameters  under  static  loading  tests.  Further,  it 
was  reasoned,  if  the  static  case  could  be  fully  modeled 
mathematically  the  extension  to  the  moving  load  case  would 
be  greatly  facilitated. 

The  laboratory  facilities  of  the  Ecole  Polytechnique 
Federale  de  Lausanne  (EPFL) ,  at  Ecublens  in  Switzerland  were 
kindly  made  available  to  the  authors.  The  first  step  was  to 
construct  of  physical  model  (Figure  4)  to  contain  soil  to  be 
tested  and  to  permit  the  implantation  of  a  regular  array  of 
lead  shot.  Loadings  would  be  introduced  at  various 
positions  on  the  soil  surface  through  a  rigid  plate.  Changes 
in  the  locations  of  the  shot  induced  by  the  surface  loadings 
would  be  monitored  by  remote  sensing  (X-rays) .  This  test 
regime  provided  the  experimental  data  for  the  project. 

Test  1 

The  test  model  was  a  glass-sided  container  with  internal 
dimensions,  80  cm  long  by  68  cm  deep  by  19  cm  thick.  The 
glass  sides  were  1.5  cm  thick  (Figure  4). 

Loads  were  applied  to  the  surface  of  the  material  through  a 
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15  cm  by  19  cm  rigid  plate.  A  sheet  of  sandpaper  was  glued 
to  the  base  of  the  plate  to  provide  a  "rough”,  frictional 
interface  between  the  load  and  the  soil  The  load  was 
varied  by  a  hydraulic  pump.  The  rigidity  of  the  load 
support  frame  was  sufficient  to  ensure  that  the  load  was 
applied  vertically  to  the  soil.  The  loading  plate  could  be 
moved  and  locked  into  position  so  that  adjacent  loaded  areas 
were  contiguous. 

Two  dial  gauges  were  employed  to  measure  vertical 
displacements,  one  on  each  side  (left  and  right)  of  the  load 
plate.  These  allowed  a  check  to  be  made  on  the  "tilt"  of 
the  loading  plate,  and  permitted  the  operator  to  judge 
whether  deformations  had  achieved  a  constant  value.  The 
criterion  used  to  gauge  final  settlement  was  that  used  by 
BOURDEAU(1),  i.e.,  Swiss  Standard  SN  670'327,  which  set  a 
settlement  rate  of  less  than  0.02  mm/min  as  a  practical  end 
point. 

Measurements  of  spatial  deformations  were  achieved  by 
placing  an  array  of  lead  shot  throughout  the  material  in  the 
model  and  registering  their  locations  at  given  times  during 
the  tests  by  exposing  the  model  to  X-rays,  the  positions  of 
the  lead  shot  being  recorded  on  photographic  film  (Fig.  5) . 

The  details  of  the  X-ray  emitter  are: 


Type: 

Phillips  MG150 

Tube: 

MoD  151/ Be 

Apperture: 

3  x  3  mm 

Focus : 

60  cm 

Film: 

Kodak  Industrex  CX 

Intensity: 

28  mA 

Power: 

120  kv 

Exposure: 

6-12  min 

Table  1. 
Technical  Data 


lead  shot 


photographic  fi 


Figure  4  Schematic 


view  of  test  mode  and  x-ray 
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Test  Material 


The  soil  used  in  the  tests  was  a  well-graded  silt  from  the 
Biolley-Orjulaz  quarry,  Switzerland,  from  which  the  gravel 
fraction  had  been  scalped.  All  standard  laboratory  tests  on 
the  material  confirmed  it  as  being  an  essentially 
cohesionless  when  dry  (Tables  2  -  4)(1). 


Gradation 

Consistency 

di  0 

(nun) 

<0.001 

Liquid  Limit 

32 

d30 

(mm) 

0.0045 

Plastic  Limit 

19.4 

d 

60 

(mm) 

0.02 

Plasicity  Index 

12.6 

* 

>20 

(d 

V  30 

)/diQ/d6o 

>1 

USCS  Classification 

CL 

d 

(mm) 

0.04 

s 

d 

(mm) 

0.10 

Composition: 

V 

d 

250 

Carbonate 

60% 

Quartz 

15 

- 

•  20% 

K 

S 

(kN/m3) 

27.1 

Muscovite 

20 

-  25% 

Table  2.  Test  Material  Particle  Characteristics 

(after  BOURDEAU) 


Direct  Shear 

moist,  unconsolidated 
saturated,  consolidated 

0=0  =29 

pic  res 

0=0  =  29 !  5 

pic  res 

m 

Simple  Shear 

saturated,  consolidated 

0  =  0  =  30.5 

pic  res 

C  =  0 

Triaxial 

saturated,  consolidated  (CD) 

saturated,  consolidated  (CU) 
saturated,  consolidated  (CU) 

0  ,=  <P  =33 

pic  res 

0=0  =33° 

pic  res 

0=0  =34° 

pic  res 

m 

Table  3.  Test  Material:  Shear  test  Results  (after  BOURDEAU) 


34 


Maxima 

Minima 

7^  (kN/m3)  17.85 

a 

e  1.27 

max 

n  56.10 

max 

r  (kN/m3)  11.92 

a 

e  0.52 

Bl  n 

n  34.20 

■  i  n 

Table  4-  Test  Material:  Compaction  Limits  (after  BOURDEAU) 

Test  materials  were  in  a  dry  condition  for  all  tests,  having 
been  dried  previously  in  an  oven  to  constant  weight. 
Moisture  content  determinations  were  taken  from  each  test 
and  indicated  moisture  contents  less  than  1%. 

Soil  to  be  tested  was  placed  in  the  test  model  in 
approximately  one  inch  layers.  The  material  was  rained 
gently  from  a  scoop  from  a  height  of  approximately  one  inch. 
The  surface  of  each  layer  was  carefully  leveled  to  minimize 
disturbance.  After  each  layer  was  placed,  lead  shot 
(average  diameter  2.5  mm)  were  placed  into  position  in  the 
center  plane  of  the  model  from  a  device  which  spaced  them  at 
approximately  1/2  inch  spacing,  while  ensuring  that  they  lay 
in  the  central  vertical  plane  of  the  model. 

When  the  model  was  full  (approx  40  cm  deep)  its  final 
dimensions  were  recorded  as  was  the  weight  of  soil.  From 
this,  its  bulk  density  was  computed.  Knowing  the  solid 
density  of  the  silt,  the  bulk  porosity  was  then  computed  and 
used  as  the  initial  reference  porosity. 

Total  Material  Weight 
Total  Material  Volume 
Bulk  Density,  r0 
Solid  Density, 

Bulk  Porosity,  n^ 


W 


Vvt 

100-  ^ys"yB) 

V(V1} 
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All  soil  porosities  computed  during  the  course  of  the 
testing  were  based  on  differences  from  the  initial  porosity, 
which  was  assumed  to  be  uniform  throughout  the  model. 

Deformation  Fvqjgtrati<?n 

After  loading,  the  movement  of  the  lead  shot  was  recorded  on 
photographic  film  after  exposure  to  X-rays.  In  order  to 
provide  fiduciary  reference,  two  series  of  lead  shot  (24 
shot  total)  were  glued  to  the  glass  walls  of  the  model: 
these,  being  at  a  fixed  position  for  every  subsequent 
exposure,  provided  a  reference  basis. 

The  first  exposure  of  each  model  test  was  made  with  the 
model  empty  and  with  the  photographic  film  held  vertically 
in  the  plane  of  the  shot  (i.e.,  midway  between  the  glass 
walls) .  This  initial  exposure  provided  the  coordinates  to 
which  the  reference  shot  images  would  be  reduced  for  each 
subsequent  exposure  throughout  the  remainder  of  the  tests. 

The  next  exposure  was  taken  with  the  model  filled,  but 
unloaded  i.e.,  under  virgin  conditions.  In  this  state  and 
subsequent  exposures,  the  photographic  film  was  held 
vertically  outside  the  model.  Subsequent  exposures  were 
made  while  the  material  was  under  load,  or  after  removal  of 
loading.  Geometric  rectification  of  these  exposures 
provided  the  coordinates  of  the  lead  shots  at  each  step  of 
the  experiment,  reduced  to  the  physical  plane  of  the  shot 
within  the  model. 

The  coordinates  of  the  images  of  the  lead  shot  were 
abstracted  from  the  developed  negatives  on  a  table 
digitizer.  The  coordinates  were  rectified  by  the 
eight-parameter  method<30)  (i.e.,  scale,  rotation  and 


36 


translation)  and  reduced  to  the  plane  of  the  shot  in  the 
model.  The  transformation  relationships  are  given  by: 

a^  +  a  *x'  +  a  *y'  a  +  a  *x'  +  a  *y' 
yH  =  0  1 _ ,  2  v”  =  3  4 _ 5-  _  (25) 

1  2  1  1  2  1 

where  the  eight  parameters,  a}  g  and  2,  are  found  by 
least  squares  regression  using  the  coordinates  of  the  images 
of  the  24  fixed  reference  shot  in  the  first  exposure 
(x",y"),  and  the  corresponding  24  image  coordinates  on  each 
successive  exposure  (x',y').  The  resulting  reduced 
coordinates  were  found  to  be  more  satisfactory  than  those 
produced  by  the  simpler  five-parameter  transformation 
initially  used  by  BOURDEAU(1). 

Exposure  times  of  the  photographic  plates  to  the  X-rays  were 
found  to  be  variable  as  the  density  and  definition  of  the 
negative  material  was  sensitive  to  the  maturity  of  the 
developing  and  fixative  solutions.  When  fresh,  exposures 
times  could  be  as  little  as  six  minutes,  while  an  exhausted 
developer  required  an  exposure  of  as  much  as  15  minutes  to 
provide  satisfactory  image  contrast. 

The  required  exposure  duration  was  a  limiting  factor  in 
attempting  to  capture  transient  effects.  However,  little 
could  be  done  to  improve  the  system.  No  other  available 
system  could  provide  the  internal  metric  equivalent  of  an 
X-ray  "movie”. 


Test  geri.es 

The  first  series  of  tests  (October  -  November  1987)  were 
inconclusive  for  a  number  of  reasons,  among  these  were: 

a.  The  image  of  the  reference  exposure  was  poor  and 
the  digitizing  process  of  the  lead  shot  images  was 
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almost  impossible  to  obtain  with  any  reasonable 
degree  of  accuracy. 

b.  During  the  course  of  the  tests,  the  loading  plate 
was  accidentally  dropped  onto  the  surface  of  the 
model,  seriously  disturbing  the  material,  and 
rendering  the  rest  of  the  first  series  invalid. 

While  the  first  series  of  tests  was  inconclusive 
from  a  numerical  data  gathering  point  of  view,  the 
experience  was  valuable  in  the  sense  of  providing 
operator  familiarity  and  "shakedown".  It  resulted 
in  a  number  of  procedural  changes  which  were 
implemented  in  the  succeeding  series  that  led  to 
the  successful  completion  of  the  testing  program. 

The  second  series  of  tests  (March  1988)  were  again  marred  by 
difficulties;  however,  two  successful  sets  of  test  results 
were  obtained.  The  STATIC  LOAD  series  of  tests  was 
satisfactorily  conducted,  while  the  initial  MOVING  LOAD 
series  was  aborted  due  to  disturbance  of  the  model.  A 
second  MOVING  LOAD  series  was  also  satisfactorily  completed. 

The  final  series  of  tests  (June  -  July  1988)  was  completely 
successful.  Both  STATIC  LOAD  and  MOVING  LOAD  tests  were 
conducted  and  the  results  were  recorded.  It  is  of  special 
note  that  the  results  of  the  two  STATIC  LOAD  series  (March 
and  June  -  July)  are  essentially  indistinguishable.  This 
observation  was  disappointing  insofar  as  no  further 
information  was  gained  from  the  second  test,  but  nonetheless 
encouraging  in  confirming  that  the  results  were 
reproducible . 


Static  laad  Series 


In  the  static  load  series  of  tests,  the  undisturbed  model 
was  initially  subjected  to  a  central,  symmetrically  placed 
load  of  low  magnitude  (3.8  kN/m2,  0.55  psi) .  The  load  was 
maintained  constant  until  deformation  had  stabilized.  This 
took  approximately  three  to  four  minutes.  An  X-ray  exposure 
was  then  taken  and  the  load  was  removed.  Then  another 
exposure  was  made.  The  load  was  replaced  and  increased  (7.2 
kN/m2,  1.04  psi)  and  held  constant  until  deformation  had 
stabilized:  an  exposure  taken.  The  load  was  then  removed 
and  another  exposure  was  made.  This  procedure  was  repeated 
for  each  load  increment. 

Moving  Load  Series 

In  the  two  MOVING  LOAD  series,  a  similar  regime  was  followed 
as  for  the  static  case,  except  that  after  the  deformations 
due  to  the  applied  load  was  found  to  be  constant  it  and  was 
translated  laterally  so  as  to  be  contiguous  with  the 
previous  load  location.  In  each  series  there  were  five  load 
locations. 

This  sequence  of  locations  was  intended  to  represent 
••snapshots”  in  time  for  a  laterally  moving  load.  The 
discretization  into  steps  was  dictated  by  the  physical  size 
of  the  footing  and  the  bolt  hole  locations  in  the  overhead 
beam.  The  time  required  to  achieve  constancy  of  deformation 
was  extremely  short.  Much  more  time  was  required  to  effect 
the  exposure.  The  time  required  to  conduct  the  sequence  of 
events  could  not  be  reduced  with  the  apparatus  available  to 
the  authors  at  the  EPFL. 

The  reductions  of  the  X-ray  exposures  for  both  the  STATIC 
and  MOBILE  test  series  were  labelled  as  follows,  and  are  so 
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identified  in  Appendix  A: 


A.  The  STATIC  series  are  identified  by  the  word 
"STATIC”  followed  by  the  magnitude  of  load,  e.g., 
38.1  kN/m2,  followed  by  the  parameter  shown,  as: 

STATIC  -  3.8  kN/m2  -  Horizontal  Deformation 

which  designation  describes  the  STATIC  series 
exposure  under  a  load  of  3.8  kN/m2,  and  the  figure 
shows  the  contours  of  horizontal  deformation. 

In  the  STATIC  series,  only  the  following  are  shown 
in  the  appendix: 

Horizontal  Deformation, 

Vertical  Deformation, 

Horizontal  Strain  (%) , 

Vertical  Strain  (%) . 


B.  The  MOBILE  load  series  was  identified  by  the  word 
"MOBILE"  followed  by  the  code  descriptor  given  in 
Table  5: 


Load 

Position 

16.2  kN/m2 

35.0  kN/m2 

0 

0.0 

1.0 

1 

0.1 

1.1 

2 

0.2 

1.2 

3 

0.3 

1.3 

4 

0.4 

1.4 

Table  5.  Reference  Sequence  for  MOBILE  Loadings 

Following  this  descriptor,  the  parameter  of 
interest  is  given  as  in  the  STATIC  series,  but 
with  these  additions: 

Volumetric  Strain  (%)  , 

Porosity. 


40 


Chapter  Six 
DISCUSSION  OF  RESULTS 

The  results  of  primary  concern  with  respect  to  the  present 
study  are  the  vertical  and  horizontal  deformations  within 
the  body  of  the  soil  induced  by  the  loadings  applied  to  its 
surface.  The  derived  strains  and  porosities  are  of 
secondary  interest.  The  validity  of  conclusions  drawn  from 
these  secondary  results  depends  upon  the  robustness  of  the 
primary  data. 


Time  effects  during  loading 

The  theories  of  LITWINISZYN<3)  and  B0URDEAU(1)  are  based  on 
the  transient  response  of  granular  materials  to  changes  in 
boundary  stresses.  In  the  present  study,  the  time  elapsed 
from  the  application  of  load  (to  the  model  surface)  to 
apparent  final  state  was  very  short,  generally  of  the  order 
of  one  minute.  Clearly,  any  transient  process  had  ceased  by 
the  time  the  X-ray  exposures  were  made.  Exposure  times  were 
of  the  order  of  6  -  10  minutes.  Moreover,  had  there  been 
any  significant  movements  of  the  lead  shot  they  would  have 
been  highlighted  by  the  apparent  elongation  of  their 
recorded  images.  None  were  seen  on  the  developed  films. 
Analyses  of  the  results  are  made  on  the  assumption  that  all 
recorded  movements  of  the  shot  were  "instantaneous".  This 
brings  into  question  the  principal  response  mode. 
Significant  transience  implies  a  "consolidation"  mechanism, 
however  the  evidence  of  the  present  study  demonstrated  that 
a  mechanism  more  similar  to  compaction  than  consolidation 
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may  predominate.  Whereas  the  classical  mechanism  of 
consolidation  is  based  upon  the  dissipation  of  excess  pore 
pressures,  progressively  transferring  the  load  from  the  pore 
fluid  to  the  solid  skeleton,  the  response  of  a  dry 
cohesionless  material  appears  to  be  instantaneous.  That  is, 
boundary  forces  in  the  compaction  mode  are  transferred 
almost  immediately  to  the  soil  skeleton.  The  corresponding 
deformations  as  a  consequence  of  the  relative  motions  of 
particles  forming  stable  configurations. 

Incrementally  moving  load 

whereas  the  loading  scheme  used  in  Switzerland  was  by 
necessity  incremental  in  space,  it  cannot  truly  (physically) 
represent  a  moving,  or  rolling,  load.  Nonetheless,  in  light 
of  the  strong  one-dimensional  deformation  response  observed 
in  the  tests  and  to  be  discussed  later,  the  representation 
is  believed  to  be  materially  similar  to  that  of  rolling 
vehicles.  That  is;  neglecting  momentum  effects,  on  the 
basis  of  the  results  it  would  be  expected  that  the  response 
of  the  soil  body  would  be  the  same  for  either  a  rolling  or 
incrementally  moving  load  so  long  as  the  amount  of  energy 
imparted  per  unit  volume  of  the  soil  body  is  the  same  in 
both  cases  in  magnitude  and  application.  The  reported  test 
results  relate  to  a  single  pass  of  the  loading  plate  over 
the  surface  of  the  model.  However,  in  a  repeated  load 
passage  scenario  the  cumulative  response  of  the  material  to 
a  rolling  or  an  incrementally  applied  loading  scheme  must 
become  less  and  less  significant. 

Vertical  and  horizontal  deformation 
Static  Loading 

It  may  be  seen  from  Figs.  A1  thru  A24  (Appendix  A)  that  the 
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vertical  deformations  are  much  more  pronounced  and 
definitive  than  the  corresponding  horizontal  deformations. 
Significant  horizontal  deformations  are  seen  to  be 
concentrated  mainly  at  the  edges  of  the  loaded  plate.  The 
maximum  vertical  deformation  is  approximately  an  order  of 
magnitude  greater  than  the  maximum  horizontal  deformation. 
It  is  difficult  without  direct  observation  to  determine  to 
what  extent  the  horizontal  deformations  are  influenced  by 
soil-wall  friction,  however  it  is  reasonable  to  expect  that 
significant  wall  friction  would  more  likely  influence  or 
even  mask  the  relatively  small  deformations  than  the 
vertical  deformation.  Consequently,  care  should  be 

exercised  for  making  any  but  qualitative  assessments  based 
on  the  reported  horizontal  deformations. 

The  "shapes"  of  vertical  deformation  curves  observed  in  the 
model,  are  not  markedly  different  from  those  predicted  by 
other  methods.  Figures  7  and  8  show  observed  and  computed 
normalized  deformations  under  the  center  of  a  single  strip 

(31 ) 

load.  The  theoretical  values  follow  from  Egorov  (Figure 
7)  and  Kandaurov  (Figure  8) .  Kandaurov  derived  a 

"granular"  stress  distribution,  similar  to  that  of  Harr,  and 
postulated  a  modified  "Hookean"  deformation,  whereas 
Egorov's  solution  was  obtained  from  pure  classical  linear 
elastic  continuum  mechanics. 

Kandaurov 's  solution  considered  a  flexible  loading  on  a 
finite  layer  of  soil  of  thickness  h:  his  solution  wass 
essentially  one-dimensional.  He  began  by  assuming  the 
constitutive  relationship  ; 

c  =  i*cr 

z  £  z 

and  thus: 

z  =  h 

A  =  4*  I  a-  •  dz 
z  E  J  z 

z  =  0 


(27) 

(28) 


Fig  7.  Comparison  Egorov  vs  Experiment 


Depth  (cm) 


Fig  8.  Comparison  Kandaurov/Experiment 
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which  under  the  axis  of  load  becomes: 

h 


where  A^  is  the  vertical  deformation  at  depth  z, 

p  is  the  magnitude  of  applied  loading, 

E  is  Young '  Modulus  of  the  soil, 

b  is  the  half-width  of  the  strip  load 

v  is  the  Coefficient  of  Lateral  Stress 

2 

X  (”-) 

and  $(x)  -  — — •  fe  2  dz  =  2*erf(— ) 

VJn  J  V2 

o 

Figure  7  shows  the  observed  normalized  deflections  from  the 
model  tests  under  the  central  axis  of  load  and  the 
corresponding  computed  deflections  using  equation  29,  for 
different  values  of  the  coefficient  of  lateral  stress,  v. 
It  is  seen  that  the  equivalent  parameter  v  is  not  constant 
throughout  the  depth  of  the  material  under  the  load.  Golden 
postulated  that,  v  is  a  power  function  of  depth,  or  v  = 
f(zn);  for  reasons  of  dimensional  compatibility  Harr 
initially  took  n  =  0  which  rendered  it  a  constant,  in  a 
later  development  (Ref.  2  Chapter  9)  he  showed  it  varying 
with  depth,  much  as  is  reported  here. 

Estimates  were  made  of  the  nu-parameter ,  v,  using  measured 
deformations  under  the  center  of  load,  A  ,  for  the  initial 

Z 

loading  of  both  MOBILE  series,  (i.e.,  16.2  kN/m2  and  35.0 

kN/m2),  and  Kandaurov's  equation  29.  These  estimates  are 
shown  in  Figure  9. 

Egorov's  solution  (in  Figure  7)  which  was  based  on  a  rigid 
strip  loading  on  a  semi-infinite  body,  is  shown  for  the 
bounding  values  of  Poisson's  ratios  of  y  =  0  and  u  =  0.5. 
It  is  seen  that  the  observed  deformations  generally  fall 


1&2  kN/m2 


Ref.  34 
X 

33.0  VN/m2 
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outside  the  prescribed  band.  However,  in  general  the  shape 
of  the  curves  are  similar  in  form  except  near  the  surface 
where  they  exhibit  opposing  curvatures. 


Load 

Deformal 

:ion  (cm) 

kN/m2 

Loaded 

Unloaded 

0.0 

0.0 

0.0 

3.8 

5.6 

5.4 

7.2 

6.0 

5.9 

11.0 

7.3 

7.3 

14.8 

8.6 

8.5 

22.0 

11.1 

11.0 

38.1 

15.8 

15.8 

Table  6.  Contact  Load  versus  Settlement 

The  reported  "growth”  of  the  bulb  of  vertical  deformation 
with  increasing  load  is  smooth  and  well  behaved,  penetrating 
deeper  into  the  model  with  each  successively  heavier  load. 
An  accompanying  increase  in  the  concentration  of  horizontal 
deformation  under  the  edges  of  load  is  noted.  In  all  cases 
no  significant  distortion  was  observed  in  the  deformation 
field,  and  symmetry  is  remarkably  well  preserved  throughout 
the  loading  and  unloading  processes.  These  suggest  that  the 
precautions  taken  to  ensure  uniform  initial  conditions  were 
satisfactory. 


Mobile  Loading 

In  the  two  MOBILE  loading  series,  the  same  general  pattern 
of  vertical  deformation  was  observed  as  in  the  single  STATIC 
tests,  but  they  were  " stretched "  horizontally.  This 
suggests  that  the  resulting  pattern  of  deformation  could 
have  been  obtained  by  using  a  single  loading  plate  taken  to 
be  as  wide  as  the  entire  (mobile)  loaded  area.  The 
distributions  of  horizontal  deformations  appear  to  be  the 
same  for  the  single  loads, as  for  the  "wider"  loading  plate 
except  that  the  concentration  of  horizontal  deformation 
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under  the  common  edge  (where  the  previous  load  would  abut 
the  current  load  location)  is  apparentlyerased:  they  appear 
to  cancel  each  other  out  and  the  soil  has  no  "memory"  of 
previous  loading 

The  peak  magnitudes  of  the  vertical  and  horizontal 
deformations  under  the  single  and  "stretched"  mobile  loading 
arrangements  are  not  materially  different.  The  developing 
pattern  of  vertical  deformation  as  the  loading  plate 
"expands"  is  clearly  one-dimensional.  This  observation  is 
in  general  agreement  with  the  models  of  FREDLUND  and 

(32)  . 

RAHARDJO  ,  who  used  the  concept  and  mechanics  of 
consolidation  to  account  for  the  deformation  characteristics 
of  granular  materials. 

The  developed  one-dimensional  pattern  of  sc «1  response  to 
the  moving  load  (Figures  A25  thru  A88)  indicates  that 
predictions  of  vertical  deformations,  may  be  obtained  by 
scaling  single  load  responses. 

Surface  deformation  ( settlements 

The  recorded  settlements  of  the  soil  surface  under  the 
loading  plate  were  plotted  against  the  magnitude  of  applied 
load  in  Figure  10.  The  plot  demonstrates  general  linearity 
over  the  full  test  range.  Note,  however,  that  the  graph 
does  not  pass  through  the  origin  (zero  load)  .  This  is 
indicative  of  an  initial  pre-loading,  or  even  a 
pre-compaction  of  the  test  material  under  its  own  weight 
prior  to  the  first  load  application.  The  overall  linearity 
of  the  graph  suggests  that  surface  settlement  may  be 
predicted  from  the  methodology  of  a  standard  "Plate  Bearing 
Test".  Using  the  observed  data  in  the  conventional  manner 
yields  an  equivalent  uncorrected  Modulus  of  Subgrade 
Reaction,  k,  of  12  psi/in. 
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(33 ) 

POULOS  gave  the  expression  for  the  settlement  of  the 
edge  of  a  flexible  strip  loading  on  a  layer  of  finite 
thickness  as: 

»,  *  SrI.l«,'h'b>  (26) 

where  p  is  the  intensity  of  loading  on  the  strip, 
h  is  the  thickness  of  the  loaded  strip, 

E  is  Young's  modulus  of  the  soil, 

4  is  Poisson's  Ratio  of  the  soil, 

b  is  the  width  of  the  strip  load, 

I  is  an  influence  factor. 

St 

He  based  this  relationship  upon  classical  arguments  which 
require  that  Young's  Modulus  be  constant,  or  at  least 
independent  of  the  magnitude  of  load.  Equation  26  was  used 
with  the  results  of  the  model  tests  to  calculate  equivalent 
values  of  Young's  Modulus.  The  results  are  shown  in  Figure 
11  for  two  values  of  Poisson's  ratio.  It  is  seen  that 
equivalent  moduli  do  vary  with  loading  and  are  not  constant. 

In  Chapter  9  of  Reference  2,  Harr  showed  that  Schmertmann's 
method  of  settlement  analysis  ,  which  was  based  on  the 

(35) 

theory  of  elasticity  and  the  results  of  Eggestad  ,  could 
be  improved  by  accounting  for  changes  in  the  material  state 
with  depth.  He  demonstrated  that  superposing  curves  of 
normalized  vertical  strain,  e^E/Ap,  for  varying  v  values,  he 
expanded  upon  Schmertmann's  single  triangle.  Strains  were 
obtained  by  applying  a  Hookean  constitutive  law  to  his 
probabilistic  stress.  It  is  quite  remarkable  that  the 
values  of  v  taken  from  the  trace  of  Schmertmann's  triangle 
provided  such  good  agreement  with  those  obtained  from 
Kandaurov's  equation.  The  overall  shapes  shown  are  in 
excellent  agreement. 
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In  Figure  12,  the  observed  values  of  v  from  the  model  test 
are  superimposed  on  a  Schmertmann/Harr  plot  from  Ref.  20. 
Once  again  it  is  seen  thatthe  overall  agreement  is  quite 
graphic.  It  clearly  demonstrates  that  the  material  state  is 
not  uniform  once  loaded  (i.e.,  v  *  const.).  In  addition, 
directly  under  the  load,  v  tends  to  be  increasingly  passive, 
while  almost  at  an  "at  rest"  (K  )  condition  is  achieved  at 
mid-depth.  This  is  consistent  with  Harr's  earlier 
methodology. 

It  is  was  observed  through  the  entire  progress  of  the 
testing  program  that  the  "rebound"  of  the  material  after 
removal  of  the  loaded  plate  is  negligible.  Some  comparisons 
are  given  in  Table  6.  The  elastic  (or  recoverable) 
component  of  the  material  response  is  hardly  differentiable. 

Of  historical  note,  plate  bearing  tests  are  used  to  provide 
predictions  of  soil  settlement  under  load,  but  they  do  not 
provide  any  indication  of  response  internal  to  the  soil 
body.  More  sophisticated  methods  use  the  results  from  cone 
or  screw  penetrometers  to  provide  profiles  of  material 
properties  with  depth.  These  more  detailed  data  may  be  used 
in  design  methods  (Terzaghi  and  Peck(36> ,  Peck,  Hanson  and 

(37)  (38)  • 

Thornburn  ,  De  Beer  and  Schmertmann)  to  provide 
engineered  estimates  of  settlement  on  and  within  the  soil 
body.  Whereas  the  Terzaghi  and  Peck  and  Peck,  Hanson  and 
Thornburn  methods  are  based  on  bearing  capacity  methodology, 
and  the  De  Beer  method  on  consolidation  mechanism,  the 
Schmertmann  method  is  based  on  the  distribution  of  strain 
under  the  load.  Harr  recognized  the  applicability  of 
probabilistic  granular  strain  rather  than  the  classical 
continuum  strain  to  Schmertmann' s  method. 
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Strains  and  Porosity 

The  strain  information  in  Figures  A1  thru  A88  were  obtained 
using  "Green's  measure  of  strain"  to  compute  the  strains 
from  the  observed  deformations.  The  "tear  drop"  pattern  of 
horizontal  strain  is  clear  in  both  the  STATIC  and  MOBILE 
cases,  and  is  appears  to  be  repeated  with  little 
modification  with  each  successive  load.  The  magnitudes  of 
the  horizontal  strains  are  small,  of  the  order  of  2%,  and  of 
limited  extent.  The  pattern  of  vertical  strains  emphasizes 
the  one-dimensional  aspect  of  vertical  deformations  reported 
earlier.  Also  strains  are  concentrated  at  the  surface  of 
the  material  under  the  load. 

The  analyses  conducted  of  deformations  and  strains  in  the 
test  model  indicate  that  granular  materials  (as  represented 
by  the  dry  silt  used)  possess  a  "memory”  that  deformation 
"persists"  after  the  removal  of  load.  This  is  evidenced  not 
only  for  applications  of  increasing  load  magnitude  at  the 
same  location  (STATIC)  ,  but  in  the  case  of  the  application 
of  a  constant  load  the  model  surface  (MOBILE) . 


-5  -4  -3  -2  -1  0 


Depth/  half  -  width ,  z/(  B/ 2) 


Figure  12  Schmertmann ' s  distribution  and  superimposed 
NU-plots  compared  with  experimental  findings 


54 


REFERENCES 


1  BOURDEAU,  P.L. ,  "Analyse  Probabiliste  des  Tassements 
d'un  Massif  de  Sol  Granulaire",  Thesis  No.  628(1986), 
Swiss  Federal  Institute  of  Technology,  Lausanne, 
Switzerland  (in  French) 

2  HARR,  M.E.,  "Mechanics  of  Particulate  Media  -  A 
Probabilistic  Approach",  McGraw-Hill,  NY,  1977 

3  LITWINISZYN,  J. ,  "The  model  of  a  random  walk  of 
particles  adapted  to  researches  on  problems  of 
mechanics  of  loose  media".  Bull.  Acad.  Polonaise  des 
Sciences  Techniques,  Vol.  11,  No.  10,  1963 

4  ROWE,  P.W.,  "The  stress-dilatancy  relation  for 
static-equilibrium  of  particles  in  contact",  Proc. 
Royal  Society,  London,  Series  A,  Vol.  269,  1962 

5  NEMAT-NASSER,  S.,  "On  Behavior  of  Granular  Materials  in 
Simple  Shear",  Soils  and  Foundations,  Japanese  Society 
of  Soil  Mechanics  and  Foundation  Engineering,  Vol.  20, 
No.  3,  1980 

6  KITAMURA,  R. ,  "Analysis  of  Deformation  Mechanism  of 

Particulate  Material  at  Particle  Scale" ,  Soils  and 
Foundations,  Japanese  Society  of  Soil  Mechanics  and 
Foundation  Engineering,  Vol.  21,  No.  2,  1981 

7  ODA,  M. ,  "A  Mechanical  and  Statistical  Model  of 

Granular  Material",  Soils  and  Foundations,  Japanese 
Society  of  Soil  Mechanics  and  Foundation  Engineering, 
Vol.  14,  No.  1,.  1974 

8  BOUSSINESQ,  J.,  "Application  des  Potentiels  a  l'Etude 
de  l'Equilibre  et  Mouvement  des  Solides  Elastiques", 
Gauthier-Villard,  Paris,  1885 

9  FOPPL,  A.,  and  FOPPL,  L. ,  "Drang  und  Zwang" ,  R. 

Oldenburg,  Munich  and  Berlin,  1897 

10  K5GLER,  D.P.,  and  SCHEIDIG,  A.,  "Baugrund  und  Bauwerk", 
Ernst  und  Sohn,  Berlin,  1938 

11  TURNBULL,  W.J. ,  MAXWELL,  A. A.,  and  AHLVIN,  R.G., 

"Stresses  and  Deflections  in  Homogeneous  Soil  Masses", 
Proc.  Fifth  Int.  Conf.  Soil  Mechanics  and  Foundations 
Engineering,  Paris,  1961 


55 


12  GARRARD,  C.M. ,  and  MULHOLLAND,  P.,  "Stress,  Strain  and 

Displacement  in  Cross-Anisotropic  and  Two  Layer 
Isotropic  Elastic  Systems",  Proc.  Third  Conf. 

Australian  Road  Research  Board,  1966 

13  Van  CAUWELAERT,  F.  "L'£lasticitd  Anisotrope  Appliquee  a 
la  Mechanique  des  Milieux  Granulaires  et  des  Roches", 
Thesis,  Swiss  Federal  Institute  of  Technology, 
Lausanne,  Switzerland,  1983 

14  DUNCAN,  J.M.  and  CHANG,  C.Y.,  "Nonlinear  Analysis  of 
Stress  and  Strain  in  Soils",  Jnl.  Soil  Mechanics  and 
Foundations  Division,  ASCE,  Vol.  96,  No.  SM5,  1970 

15  JANBU,  N. ,  "Soil  Compressibility  as  Determined  by 
Oedometer  and  Triaxila  Tests",  Proc.  European  Conf. 
Soil  Mechnics  and  Foundation  Engineering,  Wiesbaden, 
Vol.  1,  1963 

16  LADE,  P.V.,  "The  Stress-Strain  and  Strength 
Characteristics  of  Cohesionless  Soils",  Thesis, 
University  of  California,  Berkeley,  1972 

17  POKROVSKI,  G.I.,  "Investigation  on  the  Physics  of 
Soils",  Vodgeo,  ONTI,  1937  (in  Russian) 

18  KANDAUROV,  I. I.,  "Theory  of  Discrete  Distribution  of 
Stress  and  Deformation",  Izd.  VATT.  1959  (in  Russian) 

19  MULLER,  R.A.,  "Statistical  Theory  of  Stress 
Distribution  in  a  Granular  Soil  Foundation  Bed",  Bases, 
Foundations  and  Soil  Mechanics,  Vol.  4,  No.  4,  1962  (in 
Russian) 

20  SERGEEV,  I.T.,  "The  Application  of  Probability- 
Processes  Equations  to  the  Theory  of  Stress 
Distribution  in  Non-Cohesive  Soil  Foundation  Beds", 
Soil  Mechanics  and  Foundation  Engineering,  No.  2,  1969 

21  SMOLTCZYK,  H.U. ,  "Stress  Computation  in  Soil  Media", 
Jnl.  Soil  Mechanics  and  Foundations  Division,  ASCE,  Vol 
93,  No.  SM2,  1867 

22  ENDLEY,  S.N.  and  PEYROT,  A.H.,  "Load  Distribution  in 
Granular  Media",  Jnl.  Soil  Mechanics  and  Foundations 
Division,  ASCE,  Vol.  103,  No.  EMI,  1977 

23  GOLDEN,  J.M.,  "Stochastic  Models  of  Granular  Media", 
Jnl.  Engineering  Mechanics  Division,  ASCE,  Vol.  110, 
No.  EMI,  1984 


24  REIMBERT,  M.L.  and  REIMBERT,  A.M. ,  "Retaining  Walls", 
Trans.  Technology  Publishing,  Bay  Village,  Ohio,  1974 


56 


25  MARSAL,  R.J. ,  "A  Statistical  Model  for  Granular 
Materials",  Proc.  7th  ICSFME,  Specialty  Session  13, 
Mexico,  1969 

26  SWEET,  A.L. ,  and  BOGDANOFF,  J.L. ,  "Stochastic  Model  for 
Predicting  Subsidence",  Jnl.  Engineering  Mechanics 
Division,  ASCE,  Vol.  91,  No.  EM2,  1965 

27  MARSAL,  A. J. ,  "Particle  Breakage  in  Coarse  Granular 
Soils",  Proc.  7th  ICSFME,  Mexico,  1969 

28  CRC,  Handbook  of  Physics  and  Chemistry 

29  Rourke,  R. ,  "Stress  Handbook" 

30  MIKHAIL,  E.M.,  and  ACKERMANN,  F.E.,  "Observations  and 
Least-Squares",  IEP,  1976 

31  EGOROV,  K.E.,  "Distribution  of  Stresses  in  Base  under 
Rigid  Strip  Footing",  Natichn.  Issled,  Stanstiya 
Fundament stroy a ,  No.  9,  1940 

32  FREDLUND,  D.G.,  and  RAHARDJO,  H. ,  "Unsaturated  Soil 
Consolidation  THeory  and  Laboratory  Experimental  Data", 
Consolidation  of  Soils:  Testing  and  Evaluation,  ASTM 
STP  892,  1986 

33  POULOS,  H.G.,  and  DAVIS,  E.H.,  "Elastic  Solutions  for 
Soil  and  Rock  Mechanics",  John  Wiley  and  Sons,  NY,  1974 

34  SCHMERTMANN,  J.H.,  "Static  Cone  to  Compute  Settlement 
Over  Sand",  Jnl.,  Soil  Mechanics  and  Foundations 
Division,  ASCE,  Vol.  96,  SM3 ,  1970 

35  EGGESTAD,  A.,  "Deformation  Measurements  Below  a  Model 
Footing  on  the  Surface  of  Dry  Sand",  Proc.  European 
Conf.  Soil  Mechanics  and  Foundation  engineering,  Vol. 
1,  Wiesbaden,  1963 

36  TERZAGHI,  K. ,  and  PECK,  R.B.,  "Soil  Mechanics  in 
Engineering  Practice",  2nd  Ed.  John  wiley  and  Sons,  NY, 
1967 

37  PECK,  R.B.,  HANSON,  W.E.,  and  THORNBURN,  T.H.  , 
"Foundation  Engineering",  2nd  Ed.,  John  Wiley  and  Sons, 
NY,  1974 


De  BEER,  E.E.,  "Bearing  Capacity  and  Settlement  of 
Shallow  Foundations  on  Sand",  Proc.  Symposium  on 
Bearing  Capacity  and  Settlement  of  Foundations,  Duke 
University,  1965 


38 


57 

39  FENNER,  R.T. ,  "Engineering  Elasticity",  John  Wiley  and 
Sons,  NY 

40  SCHULrZE,  E.,  "Frequency  Distr  butions  and  Correlations 
of  Soil  Properties",  Statistics  and  Probability  in 
Civil  Engineering,  Hong  Kong  University  Press, 
Distributed  by  Oxford  University  Press,  1972 

41  GRIVAS,  D.A.,  and  HARR,  M.E.,  "Particle  Contacts  in 
discrete  Materials",  Abstracts  of  the  Indiana  Academy 
of  Sciences,  Engineering  Division,  1975 

42  ROHTAGI ,  V.K.,  "An  Introduction  to  Probability  Theory 
and  Mathematical  Statistics",  John  Wiley  and  Sons,  NY, 
1976 


APPENDIX  A 


A1  thru  A24  STATIC  Series 

A25  thru  A54  MOBILE  Series 

Load  16.2  kN/m2 

A56  thru  A85  MOBILE  Series 

Load  35.0  kN/m2 


STATIC  -  22.0  kN/m2  -  Vertical  Strain  (*) 


STATIC  -  38.1  kN/m2  -  Horizontal  Deformation 

-40  -35  -30  -25  -20  -15  -10  -5  0  5  10  15  20  25  30 


STATIC  -  38.1  kN/m2  -  Vertical  Deformation 

-40  -35  -30  —25  -20  -15  -10  -5  0  5  10  15  20  25  3( 


STATIC  -  38.1  kN/m2  -  Horizontal  Strain  (*) 

-40  -35-30  -25  -20  -15-10  -5  0  5  10  15  20  25  3C 


STATIC  -  38.1  kN/m2  -  Vertical  Strain  (*) 

-40  -35  -30  -25  -20  -15  -10  -5  0  5  10  15  20  25 


MOBILE  0.1  -  Horizontal  Deformation 

-35  -30  -25  -20  -15  -10  -5  0  5  10  15  20 


MOBILE  0.2  -  Horizontal  Deformation 

35  -30  -25  -20  -15  -10  -5  0  5  10  15  20 


MOBILE  0.3  -  Horizontal  Deformation 

35  —30  -25  -20  -15  -10  -5  0  5  10  15  20 


A28 


MOBILE  0.4  -  Horizontal  Deformation 


o  «r>  o  «o  o  to  o 
o  i  i  l  I  l  I  l  i 


MOBILE  0.1  -  Vertical  Deformation 

—35  —30  -25  -20  —15  —10  —5  0  5  10  15  20 


MOBILE  0.4  -  Vertical  Deformation 

-35-30  -25  -20  -15  -10  -5  0  5  10  15  20 


I 


MOBILE  0.2  -  Horizontal  Strain  (fc) 

5  -30  -25  -20  —15  -10  -5  0  5  10  15  20 


MOBILE  0.4  -  Horizontal  Strain  (*) 

5-30  -25  -20  -15  -10  -5  0  5  10  15  20 


A40 


MOBILE  0.1  -  Vertical  Strain  (*) 

-35  -30  -25  -20  -15  -10  -5  0  5  10  15  2 


MOBILE  0.3  -  Vertical  Strain  (ss) 

-35  -30  -25  -20  -15  -10  -5  0  5  10  15  2 


MOBILE  0.4-  -  Vertical  Strain  (*) 

-35  -30  -25  -20  -15  -10  -5  0  5  10  15  2 


A4  6 


MOBILE  0.3  -  Volumetric  Strain  (%) 

5-30  -25  -20  -15  -10  -5  0  5  10  15  20 


MOBILE  0.4  -  Volumetric  Strain  ($) 

5-30  -25  -20  -15  -10  -5  0  5  10  15  20 


MOBILE  0.3  -  Porosity 
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A68 


A77 


MOBILE  1.2  -  Volumetric  Strain  ($) 
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